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Abstract 

The huge spin precession frequency observed in recent experiments with spin-polarized beams 
of hot electrons shot through magnetized films is interpreted as being caused by Zeeman cou- 
pling of the electron spins to the so-called Weiss exchange field in the film. A "Stern-Gerlach 
experiment" for electrons moving through an inhomogeneous exchange field is proposed. The 
microscopic origin of exchange interactions and of large mean exchange fields, leading to different 
types of magnetic order, is elucidated. A microscopic derivation of the equations of motion of 
the Weiss exchange field is presented. Novel proofs of the existence of phase transitions in quan- 
tum Al"-models and antiferromagnets, based on an analysis of the statistical distribution of the 
exchange field, are outlined. 

1 Introduction 

Effects of ferromagnetism have been known since antiquity. But a mathematical understanding of 
the microscopic origin of ferromagnetism has remained somewhat elusive, until today! Pauli para- 
magnetism, ferro-, ferri- and antiferromagnetism are quantum phenomena connected to the spin of 
electrons and to Pauli's exclusion principle. The theory of paramagnetism in (free) electron gases 
is quite straightforward, 1 . Antiferromagnetism is relatively well understood: A mechanism for 
the generation of antiferromagnetic exchange interactions has been proposed by Anderson 2 , who 
discovered a close relationship between the half-filled Hubbard model and the Heisenberg antifer- 
romagnet using perturbative methods; (see also [3] for mathematically more compelling and more 
general variants of Anderson's key observation). It has been proven rigorously by Dyson, Lieb 
and Simon [I], using the method of infrared bounds previously discovered in [5], that the quantum 
Heisenberg antiferromagnet with nearest-neighbour exchange couplings exhibits a phase transition 
accompanied by spontaneous symmetry breaking and the emergence of gapless spin waves, as the 
temperature is lowered, in three or more dimensions. (The Mermin- Wagner theorem says that, in 
(one and) two dimensions, continuous symmetries cannot be broken spontaneously in models with 
short-range interactions, |Hj-) 

Our mathematical understanding of ferromagnetism is far less advanced. Some kind of heuristic 
theory of ferromagnetism emerged, long ago, in the classic works of Heisenberg, Bloch, Stoner, Dyson, 
Landau and Lifshitz, and others Various insights have been gained on the basis of some form 
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of mean-field theory, with small fluctuations around mean-field theory taken into account within a 
linear approximation. This approximation, however, is known to break down in the vicinity of the 
critical point of a ferromagnetic material, where nonlinear fluctuations play a crucial role, [SJ. 

In a variety of tight-binding models of itinerant electrons, ferromagnetic order has been exhibited 
in the ground state (i.e., at zero temperature); see [5|.|1U|.[TT]. One of these models is a fairly nat- 
ural two-band model in which ferromagnetism arises from a competition between electron hopping, 
Coulomb repulsion and an on-site Hund's rule, (Hund's rule says that the spin-tripled state of 
two electrons occupying the same site is energetically favoured over the spin-singlet state. It should 
be emphasized, however, that a mathematically rigorous derivation of Hund's rule in atomic physics 
from first principles has not been accomplished, so far.) None of the results in |9"].|1L)|.[TT] comes close 
to providing some understanding of ferromagnetic order and of an order-disorder phase transition 
at positive temperature. It is not known how to derive, with mathematical precision, an effective 
Hamiltonian with explicit ferromagnetic exchange couplings from the microscopic Schrodinger equa- 
tion, or a tight-binding approximation thereof, of ferromagnetic materials. But even if we resort to 
a phenomenological description of such materials in terms of models where ferromagnetic exchange 
couplings have been put in by hand we face the problem that we are unable to exhibit ferromagnetic 
order at low enough temperature and to establish an order-disorder phase transition in three or more 
dimensions. No mathematically rigorous proof of the phase transition in the quantum Heisenberg 
ferromagnet is known, to date! (Such a result has, however, been established for classical Heisenberg 
models in 0.) 

Ab-initio quantum Monte Carlo simulations of models of quantum ferromagnets are plagued by 
the well known "sign- (or complex-phase) problem". 

Thus, until now, there are neither substantial mathematically rigorous results on, nor are there 
reliable ab-initio numerical simulations of, realistic models of ferromagnetic metals, such as Ni, Co 
or Fe! Given that ferromagnetism is among the most striking macroscopic manifestations - apparent, 
e.g., in the needle of a compass - of the quantum-mechanical nature of matter, this is clearly a desolate 
state of affairs. 

In the present paper we shall not remedy this unsatisfactory situation. However, first, we attempt 
to draw renewed attention to it, and, second, we outline a formalism and some fairly elementary 
analytical observations of which we hope that they will ultimately lead to a better, mathematically 
rather precise understanding of ferromagnetism. Were it not known already, our analysis and the one 
in would make clear that ferromagnetism is a non-perturbative phenomenon involving strong 
correlations and gapless modes. To understand it mathematically will most probably necessitate a 
full-fledged multi-scale (renormalization group) analysis. The formalism presented in this paper and 
our calculations are intended to provide a convenient starting point for such an analysis. 

Analytical work on ferromagnetism may seem to be rather unfashionable. However, there are 
recent developments, such as spintronics, fast magnetic devices, etc. that may make work like ours 
appear worthwhile. Our own motivation for the work that led to this paper actually originated in 
studying recent experiments with beams of spin-polarized, hot electrons shot through ferromagnet- 
ically ordered films consisting of Ni, Co or Fe that were carried out in the group of H.C. Siegmann 
at ETH; see |12j.|13|. Back in 1998, it became clear to one of us that the concept of the "Weiss 
exchange field' (see [2], [IS]) would play a useful role in a theoretical interpretation of the experi- 
mental results reported in |12|.|13|. More generally, the Weiss exchange field appears to offer a key 
to a systematic study of phase transitions in magnetic materials, magnetic order, spin precession 
and magnon dynamics. In this paper, we focus on elucidating the microscopic origin of the Weiss 
exchange field, the role it plays in the theory of magnetism, and its dynamics. 

Our paper is organized as follows. 

In Section 2, we describe the experiments reported in [T^], an d sketch a phenomenological 
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interpretation, based on scattering theory, of the results found in these experiments, merely adding 
some conceptual remarks to the discussion of our experimental colleagues and describing the role 
played by the Weiss exchange field. We also propose some further experiments, in particular a 
Stern-Gerlach experiment for electrons traversing an inhomogeneous exchange field. 

A mathematically precise analysis of the scattering of electrons (or neutrons, photons, ...) at 
dynamical targets, such as magnetic films, metallic solids, liquid droplets, will appear elsewhere; 
(some first results appear in jTB]). 

In Section 3, we reformulate one-band t — J and Hubbard models in terms of a dynamical Weiss 
exchange field. For this purpose, the standard imaginary-time functional integral formalism for 
the analysis of thermal equilibrium states of quantum many-body systems is recalled. The Weiss 
exchange field is seen to be a Lagrange-multiplier field in a Hubbard-Stratonovich transformation 
of the original functional integral that renders the action functional quadratic in the Grassmann 
variables describing the electronic degrees of freedom; see e.g. |15j.|17|. The effective field theory of 
the Weiss exchange field is obtained after integrating over those Grassmann variables. The effective 
(imaginary-time) action functional of the exchange field and identities for Green functions of spin 
operators are derived. 

In Section 4, we determine the leading terms of the effective action of the Weiss exchange field, 
W, in the approximation where fluctuations of the length of the exchange field are neglected. For 
this purpose, we derive the ferro- and antiferromagnetic mean-field equations from the exact effective 
action of the exchange field. By solving these equations we determine the most likely length, Wq, of 
the exchange field. From that point on, the length of the exchange field is frozen to be |W| = Wq. 

We then consider a one-band Hubbard model with a half- filled band and find that, in this sit- 
uation, the effective action of W is the one of a nonlinear cr-model with a minimum that favours 
Neel order. This result is found on the basis of controlled perturbative calculations and goes beyond 
linear stability analysis of the antiferromagnetic mean-field solution, (which has been presented, e.g., 
in ^Sj). It represents a functional- integral version of Anderson's basic observations 0. 

We then turn to ferromagnetically ordered mean field solutions and show that x-independent 
fluctuations are not a source of instability of such a solution. Then we consider a one-band Hubbard 
model with a weakly filled, fairly flat band. In this situation one expects that ferromagnetism prevails. 
Indeed, we find that, at low temperatures, the ferromagnetic mean-field equation has a non-trivial 
solution, and that this solution belongs to a quadratically stable critical point of the effective action 
of W . This conclusion is the result of somewhat subtle calculations involving processes close to the 
Fermi surface, which make the dominant contribution (but would lead to small-energy denominators 
in a purely perturbative analysis). Details will appear in |18j . Our calculations support the idea 
that the one-band Hubbard model with a weakly filled, fairly flat band describes coexistence of 
metallic behaviour with ferromagnetic order, at sufficiently low temperatures. A similar conclusion 
was reached, tentatively, in for some two-band Hund-Hubbard models. The methods of the 
present paper also apply to the model discussed in (TJ; see [T%] . 

In the next to last subsection of Section 4, we exhibit a universal Wess-Zumino term in the effective 
action of the Weiss exchange field W and calculate its coefficient, which is purely imaginary. The 
Wess-Zumino term is "irrelevant" for antiferromagnets, but plays a crucial role in the dynamics of 
magnons in ferromagnetically ordered systems. Repeating arguments in JHj) we derive the Landau- 
Lifshitz equations for magnons in a ferromagnet. 

Finally, we draw attention to two well known arguments explaining why there is no magnetic 
ordering at positive temperature, in one and two dimensions; (but see jH]). 

In Section 5, we sketch novel rigorous proofs, based on analyzing the effective field theory of the 
exchange field, W, of the existence of phase transitions and magnetic order at low temperatures in 
a class of Ay-models, Heisenberg antiferromagnets and ferromagnets of localized SO(2n)- spins, for 
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n = 1,2,.... See [I] for the original results. Our proof is based on establishing reflection positivity 
of the effective field theory of the exchange field W and then using the original techniques developed 
in 0; (see also [2U], [21]). 

It should be emphasized that the concept of the Weiss exchange field has a number of further, 
quite exciting applications. We hope to return to these matters in future papers. 

Some of the material in this paper has a review character; but some of it is new. We hope 
it is fairly easy to read. If it draws renewed attention to some of the deep technical problems in 
the quantum theory of magnetism it has fulfilled its purpose. We gratefully dedicate this paper, 
belatedly, to two great colleagues and friends of the senior author (J.F.): G. Jona-Lasinio, on the 
occasion of his seventieth birthday, and H.-C. Siegmann, on the occasion of his retirement from ETH. 

Acknowledgements. J. Frohlich thanks S. Riesen and H.-C. Siegmann for very stimulating dis- 
cussions of the experiments in P. Wiegmann for some crucial advice with the calculations in 
Section 4, and the IHES for hospitality during much of the work on this paper. We all thank M. 
Azam for very useful and pleasant discussions on the uses of the exchange field. 



2 Real and gedanken experiments involving the Weiss exchange 
field 

We start this section with a brief description of recent experiments carried out by Oberli, Burger- 
meister, Riesen, Weber and Siegmann at ETH-Zurich ^21 > |13l . In these experiments, a beam of hot, 
spin-polarized electrons is shot through a thin ferromagnetic film (Ni, Co, or Fe) and the polarization 
of the outgoing beam is observed. Their experimental setup is as described in Fig. ^ 
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Figure 1: Experimental setup 

The following quantities are measurable: 

i) The thickness, d, of the film; d is a few nanometers. 

ii) The average energy, E, of an incident electron; if Ep denotes the Fermi energy of the magnetic 
film then E — Ep varies between AeV and 16eV. The group velocity of the electrons inside the 
film is denoted by v; it is not directly measurable, but it is comparable to (2(E + eV) /m* cl ) 1 / 2 , 
where eV is the average potential energy of an electron and m*j its effective mass inside the 
film. 

iii) The degree, Po, and the direction, no, of the spin polarization of the incident electron beam; 
(in Fig. ^ no is parallel to the x-direction, v to the y-direction); the same quantities, P and 
n, for the outgoing beam. 



4 



iv) The direction of the magnetization, M, of the film (in Fig. n cnosen to be parallel to the z-axis); 
the angles, Oq and 0, between no and M and between n and M, respectively (Oq = 7r/2, in 
Fig. Experimentally, the angle is found to be considerably smaller than Oq, i.e., the spins 
of the transmitted electrons rotate into the direction of the spontaneous magnetization M of 
the film. This is interpreted as being mainly due to an enhanced absorption of minority-spin 
electrons, as compared to majority-spin electrons; see (vi). (It appears that the contribution 
of spin flip processes accompanied by magnon emission into the film - "Stoner excitations" - 
to the total spin rotation is only around 5%, 

v) The spin precession angle, e, between the projections of no and of n onto the plane perpendicular 
to M (the xy-plane of Fig. [J); e is found to be "large" (tens of degrees). 

vi) Let I be the intensity of the incident beam, and let I + = I + (E) and I~ = I~(E) be the 
intensities of the outgoing beam of electrons with spin parallel or antiparallel to M, respectively, 
assuming the incident beam has intensity I and the spins of its electrons are parallel to M (no 
parallel to M, Pq = 1), or antiparallel to M (no anti-parallel to M, Pq = 1), respectively. Then 
Oq = = 0, or Oq = = n, respectively, and e = 0. I + and I~ can be measured and yield the 
spin-transmission asymmetry 



A is found to be positive and large. This is interpreted in terms of rates of transitions of 
electrons into unoccupied 3d states (holes): There are more unoccupied 3d states in the film 
with spin antiparallel to M (minority spin) than with spin parallel to M (majority spin). This 
explains qualitatively the experimental results found for A and for — Oq (see (iv)); |12j . 

vii) The orbital deflection angle, a, between the directions of the incident and the transmitted 
beam (not indicated in Fig. Experimentally a is found to be negligibly small. This tells 
us that the integrated Lorentz force on the electrons transmitted through the film is tiny. 
The precession of the spins of the electrons when they traverse the film can therefore not be 
explained by Zeeman coupling of the spins to the magnetic field inside the layer. It is mainly 
due to Zeeman coupling of the spins to what will be called the Weiss exchange field. In iron, 
the Weiss exchange field causing the observed spin precession would correspond to a magnetic 
field of roughly 8000 Tesla (which is gigantic). 

A theoretical interpretation of the experimental results reported in |12| can be attempted within 
the formalism of scattering theory. If the luminosity of the incident beam is low we can consider 
a single incoming electron. The incoming state is described as a tensor product of a Pauli spinor, 
ip- in , describing the incident electron and a state, £, of the film. Typically, £ is the ground state 
(temperature T = 0) or a thermal equilibrium state (T > 0) of the film. The outcoming state, long 
after the interactions between an outgoing electron and the film have taken place, is more complicated 
and will, in general, exhibit entanglement between the electron and the degrees of freedom of the 
film. If only measurements far away from the film are performed, as in |12j . the outgoing state can 
be described as a density matrix 



where pjy is & non- negative, trace-class operator on the Hilbert space of N outgoing electrons (the 
incident electron has knocked N — 1 electrons out of the film), N = 2, 3, . . . ; po is a non-negative 
number, the absorption probability, p\ is a non-negative trace-class operator on the Hilbert space 






(2.2) 
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of square-integrable Pauli spinors and describes the (generally mixed) state of one outgoing electron. 
"Conservation of probability" implies that 

oo 

po + Trpjv = 1. (2.3) 
JV=1 

The state P ou t is obtained by taking a partial trace of the outgoing state of the total system, including 
the film, over the degrees of freedom of the film. This is justified, because the degrees of freedom of 
the film are not observed in the experiment. If the energy, E, of the incident electron is below (or 
comparable to) the threshold, £( 2 ), for emission of two or more electrons from the film then 

p N = for N > 2. (2.4) 

In the interpretation of the experimental data provided in ^2], this is tacitly assumed. 

Experimentally, the absorption probability po and the spin polarizations Po and P of the incoming 
and the outgoing electron, respectively, are measured. The vectors Po and P are given by 

P = P n = (?/>in,<n/>in), (2.5) 

and 

P = Pn=— (2.6) 

where "0m is the wave function of the incident electron, a = (a x , a y , a z ) is the vector of Pauli matrices, 
and p\ := [Tr-^(/9i)] _1 pi is the conditional state of the outgoing electron, given that it has not been 
absorbed in the film. If Q2.4J1 is assumed to hold then 

Trn(pi) = 1 - Po (2.7) 

Since we have assumed that the incoming electron has been prepared in a pure state, 

P = |Po| = 1. (2.8) 

However, there is, a priori, no reason why p\ should be a pure state. If it were pure then 

P=|P| = 1. (2.9) 

It would be highly interesting to estimate, experimentally, the amount of entanglement with 
the film (or decoherence) in the state of the outgoing electron by measuring the quotient P/Pq. If 
P/Po < 1 then pi = [Tr^(/3i)] _1 pi is not a pure state, anymore, meaning there is entanglement with 
the film. Apparently, P/Pq has not been measured accurately, yet. 

A moment's reflection shows that spin flip processes accompanied by magnon emission in the film 
("Stoner excitations") lead to entanglement; while absorption of electrons into unoccupied 3d states 
need not be correlated with entanglement of the states of those electrons that do traverse the film. 
In fact, it is implicitly assumed in |12j that if Stoner excitations are neglected then p\ is close to a 
pure state (at least in spin-space). The experimental techniques of ^2j could be used to test this 
hypothesis. 

Next, we express the spin transmission asymmetry A (see (vi)) in terms of outgoing states. Let 
■0jt and ?/>j~ be incoming states with spin polarization Po = P^ parallel to M (majority spin) and 
Po = Pq antiparallel to M (minority spin). Let P^ ut = (p%)^ =0 and P ~ t = (p~n)n=o be outgoing 
states corresponding to ip^, respectively. Then 

a = /° ; p ° _ , (2.io) 

2 - Po ~ Po 
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as follows from ()2.3j) . More interesting would be measurements of 

P ± = Tr(pfa). (2.11) 

Clearly are parallel or anti-parallel to M; but their lengths P ± = |P ± | ought to be measured. In 
[T2"] . it is tacitly assumed that = 1, and that the states pf are close to pure states; but serious 
experimental data backing up this hypothesis appear to be lacking. It is clear that it would be 
invalidated if "Stoner excitations" played an important role. 

In the following, we outline a phenomenological description of the experiments in ^2] ; assuming 
that (|2.4|) and the hypothesis just discussed (purity of pf) are valid. (A more detailed discussion 
of the scattering approach to electron transmission- and reflection experiments will be presented 
elsewhere.) 

When an incoming electron enters the film it occupies an empty state of the film. If the film is 
crystalline this state belongs to a band of states; let a be the corresponding band index. The state 
of an electron in band a is described by a Pauli spinor 

<P=^ + A~), (2-12) 

where </> + and <f>~ are the components of (j) with spin parallel or anti-parallel to the magnetization 
M, respectively. Adopting the approximation of the Peierls substitution, the Pauli equation for 4> m 
configuration space has the form 

ihV (j) = £ a (-ihV)4>, (2.13) 
where £ a (p) is the band function of band a, 

V j = — + i%A j + i (Wj ■ - ] . (2.15) 

3 dx3 h 3 V 1 2) y ' 

Here 4> c denotes a (complex) electrostatic potential, Wfj is a (complex) Weiss exchange field, (Aj) = 
(A%, A2, A3) is the electromagnetic vector potential, and (Wj) = (Wi, W2, W3) is an STJ(2)-vector 
potential responsible for spin-orbit interactions. As discussed in |14j . (|2.13l) displays electromagnetic 
U(l)-gauge invariance and SU(2)-gauge invariance, i.e., covariance with respect to local SU(2)- 
rotations in spin space. A number of important consequences of these gauge symmetries have been 
pointed out in [Tlj . 

According to (vii) above, effects of the electromagnetic vector potential A are apparently negli- 
gible; so A is set to 0. The electrostatic vector potential (fi c is given, approximately, by 

<p c = V + iv (2.16) 

where eV is the surface exit work, and the imaginary part, v, of (fi c provides a phenomenological 
description of spin-independent inelastic absorption processes inside the film. Velocity-dependent 
spin- flip processes due to spin-orbit interactions appear to play a very minor role in the experiments 
reported in ^2]; so we may set "Wj to 0, j = 1, 2, 3. 
The Weiss exchange field W[j is given by 

= W - iw, (2.17) 

where the real part W describes exchange interactions between the incoming electron and the elec- 
tron density of the film, and the imaginary part, w, yields a phenomenological description of spin- 
dependent absorption processes. 
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Let 0i n denote the Pauli spinor describing the state of an electron when it enters the film at some 
time to- Eq. (|2.13l) can be solved for (f> = <j) t , t > to, with (j) to = (pin- The solution is explicit if A = 0, 
Wj = 0, j = 1, 2, 3. Let us suppose that the real part W and the imaginary part w of the exchange 
field Wq are both anti-parallel to the magnetization M of the film. Then ()2.13j) leads us to consider 
two simple, quasi-one-dimensional scattering problems in a complex potential well of depth 

eV + iev± -(n~iu), (2.18) 

with O = |W| and u> = |w|, for electrons with spin parallel to M(+), or anti-parallel to M(— ), 
respectively. The solution to these scattering problems can be found in every book on elementary 
quantum mechanics. For the purposes of interpreting the results in ^2], a semi-classical treatment 
appears to be adequate. The group velocity, v±, of an incoming electron wave with energy peaked 
at E and spin up (+), or down (— ), inside the film can be found by solving the equation 

£ a { V ) = E - eV ^hl (2.19) 

for p and then setting 

v± = ^(p*), (2.20) 

where is a solution of (|2.19|) chosen such that v± points in the positive y-direction. The sojourn 
time t± of the wave inside the film is then given by 

d 

t± = — , with v± = v± . (2-21) 

v± 

If 0^ is the state of the electron with spin up (+), or spin down (— ), respectively, when it enters the 
film its state 4> ovA upon leaving the film is then given, approximately, by 1 

0± t - exp {-ir ± [{E - eV) - iev T ^(O - iu>)]\ <fc. (2.22) 

The presence of inelastic absorption processes implies that ev =F w/2 > 0. Setting 

exp(-r ± [{)T^/2]) = cy/1 ± A, (2.23) 

we find that 

^ ut = cVT±Ae^ 2 Ut (2-24) 

where A is the spin transmission asymmetry (see point (iv) and (|2.1(jp ). e is the spin precession angle, 
and 6 is an (unimportant) spin-independent phase. If v + = V- = v then r + = r_ = d/v, and (|2.22f) 
yields 

e^T±fi^(d/v)|W|. (2.25) 

Thus, measuring e and d and estimating v yields an approximate value for the size of the spin 
precession angular velocity and hence of the size of the Weiss exchange field. 

Eqs. (|2.22j) and (|2.24|) could be mistaken for equations describing kaon- or neutrino oscillations 
and are analogous to the equations describing the Faraday rotation of light traversing a magnetized 
medium. 



1 Here and henceforth, we use units such that h — 1. 
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The origin of the Weiss exchange field W is hardly a mystery: The spin of an electron traversing 
the film in a band a apparently experiences (exchange) interactions with the spins of the occupied 
states of the film. Since states with spin up and with spin down are occupied asymmetrically 
(corresponding to the fact that M 7^ 0), the net spin density, S(x), at a point x in the film is 
different from zero. The Weiss exchange field W = W(z) is given by 

W(x) = -JaS(aj), (2.26) 

where J a is the strength of the exchange coupling between spins in the a th band and those in the 
occupied band. The theoretical discussion above is based on a mean field ansatz: The exchange field 
W in (gJH) , (f2~T7)) is chosen to be 

W = (W(x))=-J a (S(x)). (2.27) 

(Our conventions are such that (S(x)) is parallel to the magnetization M.) One of the surprising 
implications of the experimental results of is that, apparently, J a is quite large, even for rather 
high- lying bands (a), implying that the orbitals of states in such bands must have substantial overlap 
with those of states in the partially occupied, spin-polarized band. 

Eq. l)2.26j) niakes it clear that W(x) is a dynamical field. One of the main purposes of this 
paper is to derive the effective quantum dynamics of W(x) within a Lagrangian functional integral 
formalism and to sketch what can be accomplished with this formalism. 

Before proceeding with this program, we outline some gedanken experiments in the spirit of the 
experiments in |12| . 

i) By applying an external magnetic field to the film rotating in the xy-plane with angular velocity 
ujq, the exchange field W could be made to rotate around the z-axis: 



W(t) = W (ecos(u t + d),esin(u t + 5), \/l - e 2 ) . (2.28) 



A polarized beam of electrons shot through this film must exhibit Block spin resonance; but, 
in this experiment, it would be due to the rotation of the exchange field. 

ii) One may envisage a Stern Gerlach experiment for electrons. One would start by constructing 
a sandwich of two ferromagnetic metals, I (e.g. Fe) and II (e.g. Ni), with exchange fields W/ 
and W77 of different strength, joined by a transition region, a mixture of I and II, of width do. 
The transition region would be parallel to the xy-plane (see Fig. [2J). One shoots an unpolarized 
beam of (not very hot) electrons through the sandwich along the transition region between I 
and II, as shown in Fig. [2 One would expect to detect two beams emerging on the other side 
of the film in slightly different directions that are spin-polarized in opposite directions. 

The force, /, in the z-direction on an electron with spin up/down inside the film is given, 
approximately, by 

/ ± l —+ (2.29) 

2d K J 

It yields a change in the z-component of the momentum of the electron during its passage 
through the film given by 

Ap z * ± l 1 m -, (2.30) 
2do v 

where v is the average group velocity. The deflection angle a is found from 

tana = — — . (2.31) 
Ap y 
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Figure 2: Stern-Gerlach Experiment 

Of course, as discussed above, the intensity, I + of the upper beam can be expected to be much 
larger than the intensity, I_ , of the lower beam, due to spin-asymmetric absorption inside the 
film. In order to achieve adequate focussing, the incident beam could come from the tip of a 
scanning tunnelling microscope. 

iii) In an experimental set-up similar to the one above, one would force an electron current, spin- 
polarized in the z-direction, through the film. Then a Hall tension in the z-direction should be 
observed (Hall effect for spin currents; see [Tl]^. 

3 The Weiss exchange field in general one-band (t—J and Hubbard) 
tight-binding models 

In this section, we study the origin and dynamics of the Weiss exchange field in simple tight-binding 
models. We start by considering one-band models; but it is straightforward to include higher (par- 
tially occupied or empty) bands relevant for the study of electron transmission through magnetic 
films. 

Every site, x, of a lattice T e = Z d , d = 1,2,3,... carries a four-dimensional state space, C 4 , 
corresponding to an empty state, a one-electron state with spin up (+) or down (— ), or a state 
of two electrons in a spin-singlet state. The corresponding basis vectors of C 4 are denoted by 
{|0)x, |+)xj \~)x, | H — )x}i x S T. We introduce electron creation- and annihilation operators, c\(x) 
and c s (x), respectively, where s = ±, x € T, satisfying canonical anti-commutation relations 

{c#(x),c*(x')} = 0, {ct(x),c s/ (x')} = 5 ssl S xxl , (3.1) 

with c*=corct Then 

c s (x)\0) x =0, for s = ±,xeT (3.2) 

and 

\s) x = cl(x)\0) x , etc. (3.3) 

Number operators are given by 

n s (x) = c\(x)c s (x) , n(x) = n + (x) + n_(x) , (3-4) 

where n s {x) measures the number of electrons at site x with spin s. By (|3.1|) . n s {x) has eigenvalues 
and 1, while n(x) has eigenvalues 0, 1 and 2. The spin operator, S(x), at site x is given by 

S ( x ) = \ ^Z c l( x ) a ss'C s '{x) (3.5) 

s,s' 
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where a ss i is the (ss')-matrix element of the vector of Pauli matrices a = (a x , a y , o~ z ). The operators 
S(x) vanish on lO)^ and | H — ) x , while, for any linear combination (p{x) of 1+)^ and |— ) x , 

S(x) ■ S(x)cf>(x) = ~J>(x) . (3.6) 

The dynamics of a gas of electrons moving on the lattice T is assumed to be generated by at — J 
model Hamiltonian of the form 

H = T + U + E^, (3.7) 
where the kinetic energy operator, T (the hopping term), is given by 

r = ^t(x-y)^ct(x)c s (y), (3.8) 

x^y s 

U = U J2 n +( x ) n -( x ) ( 3 - 9 ) 

X 

describes an on-site (Hubbard) repulsion, and 

= ~\ J A* ~ y)S(x) • S(y) (3.10) 

describes (effective) exchange interactions between the spins of electrons. In (|3.8|) . t(x — y) is the 
amplitude for hopping of an electron from site y to site x. Selfadjointness of the Hamiltonian H 
implies that t(x — y) = t(y — x). The constant Uq is a measure of the strength of the on-site repulsion, 
and J^(x — y) is the exchange coupling between spins at x and y, with J^(0) = 0. 

The term U can be rewritten as follows. When x is empty or singly occupied Uon+(x)n^(x) 
vanishes; when x is doubly occupied it takes the value Uq. Thus, by (|3.6|h 

U n+(x)n-(x) = \u n{x) - ^U S(x) ■ 8(x) = ^U : 8(x) ■ S(x) : , (3.11) 

where : . : denotes normal ordering. We now set J(0) = |f/o, J(%) = J^(x), for x^O, and find 

H = T + E, (3.12) 

with 

E=^J(x-y):S(x)-S(y): . (3.13) 

%,y 

The Heisenberg model corresponds to a completely flat band, i.e., 

T = 0, or i(x) = 0. (3.14) 
Our goal is to study the grand-canonical partition function 

E((3,»):=Tr(e • (3.15) 

and the thermal equilibrium state given by the density matrix 

P^ = ~(f3,n)- l e * (3.16) 
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for the family of Hamiltonians introduced above. 

Remark. Of course, these objects must first be calculated for arbitrary bounded regions, A, 
in the lattice T (indicated by a superscript "A"). Afterwards, one may attempt to pass to the 
thermodynamic limit, A f T, setting 



and 



Pg((3,^ = - lim— lnH A (/5,/i) 



^(•)=]imTr(P^(-)) 



(3.17) 



(3.18) 



where |A| is the number of sites in A. These are standard matters (see e.g. [21]) and will not be 
discussed any further. The superscript "A" will often be suppressed in our notation. 

It is convenient to study H(/3, fi) and the state ujp ifJi by making use of functional integration. With 
every x £ T and every imaginary time t € [0,/3), we associate anticommuting (Grassmann) variables, 
C s (x,t) and C s (x,t) with 

(-) (-) , , 

{C s (x,t), C s ,(x',t')} = 0, 



and 



C#(x,t + (3) = -C#(x,t). 



(3.19) 



The anti-periodic boundary conditions in (|3.19j) are a consequence of the KMS condition; see e.g. 
|17j . We introduce an (imaginary-time) action functional 



Sj(C,C) 



(It 



£ C s (x, t)^C.(x, t) j - H(C(.,t),C(.,t)) Cs(x, t)C.(x, t) 

\ / x,s 



(3.20) 



where H(C,C) is obtained from the Hamiltonian H by replacing the operators c s (x) and ct(x) 
by C s (x,t) and C s (x,t), respectively. Because of the presence of products of C- and (7-fields at 
coinciding sites and imaginary times the quartic term needs to be regularized; (the ambiguity of 
the quadratic terms is an unimportant constant). We introduce an infinitesimal separation of the 
imaginary times, replacing S(x) in (|3.13|) by 



S € {x,t) = ^ C s (x,t + e)a ss >C s > (x, t - e) 



(3.21) 



in the action. Then the time ordering of the C- and C-fields corresponds to the normal ordering of 
the operators S. Because of ([3.13)1 the regularized action is now given by 



S'j(C,C)= f dt {y^C s {x,t)^-C s {x,t)\ -^2i{x-y)^2c s (x,t)C 8 (y,t) 

V x,s / x^y s 



Berezin's integration form for anticommuting variables is formally given by 

VCVC = U mm U x U s dC s (x,t)dC s (x,t) . 



, (3.22) 



(3.23) 
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Then 



= limy VCVCe s J {C ' C) (3.24) 
and the state up ifl can be reconstructed from the imaginary-time Green functions 

(n,-<7 s . (xj^C,,. (4, ty)f} 7lM = E((3, /i)- 1 lim f VCVCe^^UjCs. (x^Cj. (z£, t£) ; (3.25) 

see, e.g., [pi, |23]. 

If Sj(C, C) were quadratic in C and C the Berezin integrals in (|3.24j) . I|3.25j) could be evaluated 
and expressed in terms of determinants. The only contribution to Sj(C,C) not quadratic in C, C 
comes from the term 

P 

fdtJ2 J(x ~ v)S e (x, t) ■ &(y, t) . (3.26) 



o 



By introducing a Lagrange-multiplier field ("Hubbard-Stratonovich transformation") this term, 
too, can be rendered quadratic in C and C. The Lagrange multiplier field will turn out to be the 
Weiss exchange field W. 

Let K denote the matrix inverse of J, i.e., K(k) = J(k)^ 1 , where A; is a point in the first Brillouin 
zone, By, of the lattice T, and J, K denote the Fourier transforms of J and K. If Uq is chosen to be 
large enough then J and K are positive-definite, i.e., J(k) > 0, K(k) > 0, for k £ Br- We introduce 
a real, vector-valued random field, ~W(x,t), x 6 T, t £ [0,/3), with 

W(x,t + (3)=W(x,t) (3.27) 

(KMS condition for bosons; see e.g. ^Zl), assumed to have a Gaussian distribution 

c^k(W) = e SK{w) VW , (3.28) 

where 

P 

Sk(W) = ~ I dtJ2 W(x, t)K{x - y)W(y, t) , (3.29) 

<p,y 

DW = const^ U te[0>p) Il x d 3 W(x, t) , (3.30) 

and consti<- is chosen such that 

Jdfi K (W) = l. (3.31) 



We define 2 



and 



S (C,C) = S J=0 (C,C), (3.32) 

P / \ 

S^C, C; W) = - J dt \J2 W ( x > *) • S ( x > *) ) « ( 3 - 33 ) 

n \ x / 



S(C,C;W) = S (C,C) + Si(C,C;W) + S K (W) . (3.34) 



"Henceforth we suppress the regularization indicated by e. 
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We note that S(C, C;W) is quadratic in C and C and in W (separately), and that 

J DWe S ^ = e^^ c ) , (3.35) 

as follows by quadratic completion. 

It is apparent from ()3.33j) that W(x, t) plays the role of the Weiss exchange field. In order to 
calculate the grand canonical partition function H(/3, fj,), we have to integrate the R.S. of (|3.35|) over 
C and C. Since 5((7, C; W) is quadratic in C and C, it is tempting to interchange the W - and the 
C, C - integrations. First doing the C, C - integration yields 

J V cvCe s ^ +s ^ C ' C ^ = detpw) , (3.36) 

where the operator Dw acts on the space 

l) = C°°([0,P)) ap ® (l 2 (T)®£ 2 ) (3.37) 
( "ap" stands for "antiperiodic" ) and is determined by requiring that 

P 

S (C,C) + S 1 (C,C;W) = J dtY,C(x,t){D w C)(x,t). 

x,s 

Thus, _D\y is given by 

d 1 
= -Q t ® 1 + 1 ® (t - Ml) + g W ( a; ' *) • °" ' ( 3 - 38 ) 

where (th)(x) = Yl K x ~ V)h{v) > ^ G ^(r). 

After Fourier transformation in t and x, vectors in the space i) are given by two-component 
spinors 

7+(Mo) 
/-(Mo) 



<^(k,/so) 

with k G By and A; = (n/P)(2n + 1), n 6 Z, and 

(I>w^)(k,*k)) = (-*Ato + t(k)-M)^(k,*o) + i y ^3^W(k-k',fco-^)-(^)(k',^), (3.39) 

where t(k) denotes the Fourier transform of t(x). Defining the effective action, S e ff(W), of the 
exchange field, W, by 

5 e// (W) = lndet(D w ) + S K (W) , (3.40) 

we find that, for example, 

3(fi,n) = j VWe s -ff {w) . (3.41) 

Note that 

_^_ e 5 ( c,c ; w) = _ (j2K( X -y)W(y,t) + S( X ,t)\ e § ^). (3 . 42) 
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Setting 

6w ^ i) = 5w(x t) + Y. K ^ x ~ y) *) > 

y 

n 

the imaginary-time Green function corresponding to a product, Y[ S a (xj,tj), of spin operators 
turns out to be given by 

n „ n 

(n^i^wH-imM -1 / fwimf^t,.)^^, (3.43) 

i=i J j=i 

where we have used 1)3.40(1 and (|3,42|) , 

Our purpose is now to determine S e ff(W) as explicitly as possible, in order to get some insight 
into the R.S. of (I3~1T1) and (I3~i3l . 

4 Effective cr-models of the exchange field 

If we want to apply the method of steepest descent to estimate the R.S. of 1)3.41(1 and of 1)3.43(1 . we 
must look for the (absolute and local) maxima of S e ff(W). This leads us to consider mean-field 
ansatze for configurations, W, corresponding to local maxima of S e ff(W). 



Ferromagnetic mean-field theory 

We set 



W(s,t) = W n, (4.1) 



independently of x and t, where n is a unit vector, and W$ > 0. Because of the symmetry of the 
problem under rotations of W, we can choose n = n z to point into the z-direction. The problem of 
evaluating 

g m f(P,KW ) = hn^-^L^W = W n z ) (4.2) 

is elementary, because, in momentum space, all modes decouple from one another, and we can appeal 
to well known calculations (see e.g. )15j . |1 7j ) to obtain the answer: 

9™I(P,K W ) = ±K(0)W§ - lim -L lndet(Dw=w ») 



1 + e -P{<K)+W /2)\ x 

x (i + e H3(*M-Wb/2))] + const., (4.3) 



where 

e(k)=t(k)- M . (4.4) 

Formula 1)4.3)1 can also be derived directly from the definition of lndet(Z?w) = trln(Z>\y) by using 
the Poisson summation formula. 
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In the limit (5 — > oo, expression Q4.3|) converges to 



g m f(oc,^W )=±k(0)Wi 

dk 



+ 

Br 



(2vr) s 



(4.5) 



+ const . , 



which is the ground state energy density of electrons in a constant exchange field W = Won 2 . The 
value of Wq minimizing g m f is found by solving the equation (d/dWo)g m f(P, \x\ Wq) = 0, i.e., 

K(Q)Wq = ~ [7% UeP^)-^) + l)- 1 - ( e /^«+^) + I)" 1 ] . (4.6) 
2 J {2ny L J 
Br 

When j3 — > the R.S. of (|4.6|) tends to 0, and we conclude that, for small /?, this equation only has 
the trivial solution, Wo = 0. When (3 — > oo, equation (|4.6|) yields the equation 

Besides the trivial solution, this equation also has a non-trivial solution Wq > 0, provided K(0) is 
sufficiently small (depending on t). We recall that 

K(0) = J(0)- 1 = ^ + J^(0)) . (4.8) 

Thus, for f/o/t*, with i* := max|i(x)|, large enough, equation (|4.6|) has a non-trivial solution, 
provided /3 is large enough, even in the Hubbard model, where = 0. Note that, by (|3.42|) and 

- £ if (x - y)(W(y, t)) AA1 = <S(x, t)) Ptli . (4.9) 
y 

Thus, in mean-field theory, 

K(0)W =M, (4.10) 
where Wq is the non-trivial solution of equation (|4.6|) . and M is the spontaneous magnetization in 

)S) 

09 - oo), 



mean-field theory (i.e., -Mn 2 = (S(x,t))"^l). Thus, equation (|4,7j) tells us that, at zero temperature 



2 y (2tt) 3 

where J(0) = ^ </(x) = if (O)^ 1 . This equation has an obvious interpretation apparent from Fig. |21 

X 

Matters simplify further for a flat band, t = 0, i.e., for the Heisenberg model. Then eq. (j4,7j) 
reduces to 

M= fl/2, -J(0)/4</i<J(0)/4, 
\ , otherwise , 

for P — > oo. For, eq. (|4.5j) implies that 

^(oo, M ; Wb) = ^(0) W 2 - ^9 - |a*|) , (4.13) 
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Figure 3: Spin bands 

up to an unimportant constant. Thus, for \fi\ < 1/(4K"(0)) = J(0)/4, g m f(oc, fi; Wo) has a quadratic 
minimum at Wo = l/(2if (0)), with 

g m noc,^W ' 



2K(0) 



and 

d 2 g m f 




oo,^W -^—\=K(0). (4.14) 



dW $ 2 ^ (0) 



Antiferromagnetic mean-field theory 

Antiferromagnetism and Neel ordering in Hubbard models is discussed in some detail in ^S], so we 
shall be brief. To simplify matters, we consider the lattices r = Z d , d = 1,2,3,.... Our mean- field 
ansatz configuration for the exchange field W is given by 

W(x,t) = (-l)HWbn r , (4.15) 

which describes Neel order. We skip detailed calculations and give the result for the free energy 



g ™f(P l ^Wo) = \k{^)Wo 2 --p J ^3 In 



1 + e -/Ve(k) 2 +(i/4)iy 2 \ ( 1 + e f3^7WMl/W$ 



const , (4.16) 



where = (n, it, . . . , ir) denotes a corner of the Brillouin zone. In order to derive (|4.16j) . we 
assume that e(k + k„-) = — e(k), i.e., we assume the band to be half-filled. Note that small deviations 
<5e(k) = e(k+k 7r )+e(k) from half filling can be taken into account with the help of perturbation theory 
in <5e(k). The idea behind the derivation of (|4.16jl is to consider simultaneously the contributions 
of k and k + k^, diagonalising a 4 x 4 matrix. The energy eigenvalues, for a fixed k, are given by 
e±(k) = ±y / e 2 (k) + (1/4) Wq, each one with twofold degeneracy. From (|4.16|) we obtain the mean- 
field value of Wo by solving the equation dg m f(f3, /i; Wq)/8Wo = 0. In the limit where (3 — > oo, this 
equation reduces to 

K(K) = - I , , 1 • (4.17) 

K } 4 J (2vr)3 ^ e(k )2 + (i /4) ^ V > 
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The R.S. diverges logarithmically when Wq — ► 0, for a lattice T of arbitrary dimension. Thus (|4.17f) 
has a non-trivial solution Wq > for arbitrarily large values of K{k^) = (AUq/3 + J^kvr))^ 1 , i.e., 
for arbitrarily small values of the on-site repulsion Uq (even if = 0), and hence the same is true 
for finite, but sufficiently large values of (3. Thus 

|e±(k)|>~Wo>0, (4.18) 

and we conclude that a strictly positive energy gap Wq > opens at the Fermi- surface. 
We note that the solution of (|4,17|) behaves like 

Wo oc Ue- ak{ ^ )u , (4.19) 

for some constant a. For the half-filled Hubbard model (J^ = 0, Uq > 0), ^(k^) = 3/(AUq), and we 
find that Wq has an essential singularity at Uq = 0. 

Comparing these results with those obtained for the ferromagnetic mean-field ansatz, we conclude 
that, at zero temperature, in the Hubbard model with K(0) = K(k n ) = 3/(4C/o), a non-trivial 
ferromagnetic mean-field solution does not exist when Uq is small, while an antiferromagnetic mean- 
field solution exists for arbitrarily small values of Uq at half filling, where the (mean-field) ground 
state energy of the Neel state is found to be below the ground state energy of the ferromagnetic state. 

The magnetic cr-model 

As we have seen, these ansatze correspond to quadratic (local) maxima of the effective action as a 
function of the length Wq of the exchange field. This suggests that the effective action S e ff(W) is 
maximal on exchange-field configurations of constant length, |W(x,i)| ~ Wq. It may therefore be 
appropriate to study the effective action ()3.4Uj) with the constraint 

| W(x, t)\=W , for all t £ [0, (3) , iGT. (4.20) 

To study magnetic order, we consider the following ansatz for the effective action on configurations 
of exchange fields of constant length, |W(x,t)| = Wq: 

P 

5 e //(W) = /^E{-^E t)K{x - y)W(y, t) + C t \d t W(x, t)\ 2 + C v \ VW(x, t)\ 2 ] 

o x y 

+ Cwz J dtds W(x, t, s) ■ (dtW(x, t, s) A d s W(x, t, s)) + higher order terms , 

S 2 ' N x 

(4.21) 

where we have set W(x, t) := H / _1 W(x, t) and where, for each x £ T, W(x, t, s) is chosen to be a 
smooth extension of W(x,i) from the circle {t|0 < t < /?} to the northern hemisphere, S 2 ' N , of the 
2-sphere of radius /3/(2tt). The operator V denotes the finite difference gradient on the lattice, and 
the corresponding term in the effective action has to be interpreted as 

J]|VW(x,t)| 2 := |W(y,t)-W(x,t)| 2 . 
x \x-y\=\ 

The functional S* e //(W) is the (imaginary-time) action of a non-linear a-model with SO(3) 
symmetry and with a Wess-Zumino term. This model is called magnetic a-model. It is expected to 



18 



describe exchange field configurations W in thermal equilibrium at inverse temperature j3, provided 
the coefficients Cf, Cy and Cwz are chosen in such a way that 

S eff (W) « S eff (W)\ lw{Xjt)l=Wo + const . , (4.22) 

on configurations W(x, t) satisfying the constraint |W(x,t)| = Wo, with S e ff(\V) as in (|3.4UI) . More 
precisely, if W = W m j + X- where W m j is a mean (-exchange) -field configuration and X 1S a small 
perturbation of W TO /, i.e., x{ x i -L W m j(x) and \x{ x i 01 ^ 1> f° r all a; € T and all times i, then we 
want (|4.22l) to hold exactly to second order in X- Eq. (|4.21|) makes it clear that, in order to calculate 
Ct, Cy and C^yz such that Q4.22JI holds to second order in x, it suffices to consider time-independent 
exchange fields, W(x,i) = W(x), to calculate Cy, while it suffices to consider site-independent 
exchange fields, W(x,t) = W(i), when calculating C% and Cvkz- Obviously, the signs of Ct and Cy 
will determine the stable mean-field configuration W m y. For W m j to be independent of time t, Ct 
must be negative (with our sign convention). For ferromagnetism to prevail, Cy must be negative. 

In the next two subsections, the coefficients Ct and Cy are determined in such a way that (|4.22|) 
holds to second order in for an appropriate choice o/W m ^ (depending on our choice of parameters 
in the original model). The coefficient Cwz of the Wess-Zumino term will be calculated and shown 
to be imaginary in the next to last subsection of Section 4. The Wess-Zumino term is crucial in 
arriving at the correct magnon dynamics for ferromagnetically ordered magnets. 

Stability of antiferromagnetic ordering 

In this subsection we calculate the coefficients Ct and Cy of the magnetic a- model (|4.21|) , in such a 
way that (|4~22|) holds, in a regime where these coefficients can be calculated perturbatively. For the 
Hubbard model this is the regime where the band is half-filled and the signs of the coefficients Ct and 
Cy will indeed determine the stable mean exchange field to be time-independent and antiferromag- 
netically ordered. 

We first derive an innocent looking, but important identity, which will also be useful in the next 
section. We set 

D w = -d t + e-W. (4.23) 

Equation (j3,38|) shows that 

£>W = ~ d t + e + ^lffi - W 2 a 2 + W 3 a 3 ) , 

where the superscript "T" stands for transposition, (i.e., taking the adjoint, followed by complex 
conjugation). Hence 

a 2 D^a 2 = -d t + e-W = D W . (4.24) 

Furthermore, 

D^ = -d t + e + W = DZ W • (4.25) 
Since a 2 = = cr^ 1 , it follows from (|4.24|) that 

det(Dw) = det(D^) = det(a 2 D^a 2 ) = det(D w ) , (4.26) 

and hence 

det(Dw) = det(Z^v) = det(ZT w ) = det(£>_ w ) • (4.27) 
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This identity shows that the real part of the logarithm of the determinant of Z?w is given by: 

5Rlndet(£> w ) = -lndetpw-D-w) , 

W 2 

D W D_ W = (8 t + e) 2 - -± - [e, W] - (8 t W) . 
We are tempted to treat [e, W\ + (dtW) as a perturbation of 

A := (d t + e) 2 - -± . (4.28) 

This is justified if |e| <C Wq and <C Wq. In this regime the band of electrons with spin parallel 

to the background field, W m y, is filled and the band of electrons with spin anti-parallel to W m j 
is empty, i.e., the lattice is half-filled. At half-filling, we expect the Hubbard model to exhibit an 
antiferromagnetic phase at sufficiently low temperatures. We recall that lndet = trln and expand 
the logarithm up to quadratic order in [e, W]/A and (dtW)/A. Terms linear in W vanish at critical 
points of the effective action; W-independent terms are omitted. We are then left with 

Sftlndet(D w ) = --tr(A _1 ([e, W] + (W))A _1 ([e, W] + (d t W))) + h.o. , (4.29) 

where h.o. stands for terms of higher order in | e | / Vl^o or |^W]/Wa. Next, we pull A -1 through 
[e, W] + dtW using the identity 

[A" 1 , [e,W] + (d t W)] = -A" 1 [A, [e,W] + (WJA" 1 , (4.30) 

which shows that the commutator only contributes higher-order terms to (|4.29|) . Mixed terms with 
one spatial and one time derivative of W do not contribute to the trace. We thus arrive at 

Klndetpw) = -^tr(A" 2 ([e,Ty] 2 + (d t W) 2 )) + h.o. . (4.31) 
Comparison with (|4.21|) then yields an expression for Cy, 

c v -S? I dk| V * )|2 ^ £ (HW-^oW «sk / " k| V£(k)|2 ' 

B T fc ° Br 

which is positive. Thus the real part of the cr-model action (j4.21|) for the Hubbard-model exhibits 
its maximum (w.r.t. time-independent exchange fields) on configurations for which 



p 

J dt^2\VW(x,t) 



is maximal, and this is the case for a Neel ordered (staggered) exchange field. 
From (|4.31|) we can also read off the Cf-coefficient in (|4.21jl to be given by 

4 8 0^ ((_^ + e(k)) 2 -Wo 2 /4) 2 4^ ' 

which is manifestly negative. From these calculations and the mean field equation Q4.17JI we conclude 
that the Hubbard model at half filling exhibits a stable antiferromagnetic phase if the on-site repulsion 
Uq is large enough as compared to the hopping term, and if the temperature is sufficiently low. 
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Linear stability of the ferromagnetic mean exchange field 

In this subsection we analyze conditions that imply stability of ferromagnetically ordered exchange 
fields. In the Hubbard model, we expect a stable ferromagnetic phase to prevail when the band is 
weakly filled and fairly flat. In such a regime, since there is no energy gap at the Fermi surface, the 
operator A defined in (|4,28|) has zeros and the perturbative expansion of the last subsection breaks 
down. We therefore have to resort to a linear stability analysis around a ferromagnetically ordered 
mean exchange field. We will see that the stability of the ferromagnetically ordered exchange field 
is due to contributions close to the Fermi surface. 

The stability of a ferromagnetically ordered exchange field w.r.t. x-independent fluctuations can 
be shown quite easily, since the operator D\y is diagonal in momentum space when W(x, t) = W(i) 
is x-independent. We can then write 



1 f cZp { 

VT JK\ ln d6t(jDw) = J (2^3 ln det ( d t + e "(P) 

Br 



+ 



W(t) 



(4.32) 



In the operator formalism outlined at the beginning of Section |31 

W(t) 



det d t + e(p) + 



where 



9 t C7 e>w (i,s) 
E/ e ,w(M) 



const xTr C 4C/ ej w(/3, 0) , 



(4.33) 
(4.34) 



and 



h e ,w (t) 



( o 



V o 





-V Ml , W(t) 

e(p)i 2 + a 




\ 




2e(p)/ 



The trace of the propagator U €j w(/3, 0) is 

Tr C 4[/ £iW (/?,0) = 1 + e-^P)T¥ c2 *7 e= o, W (AO) + e~ 2 ^ 
Using the Holder- inequality for traces (see e.g. HJ), we find that 

I Tr c2 (C7 w(/3,0)) I < Tr c2 (U 0tW=Wo . a {/3 , 0)) , 



(4.35) 



(4.36) 



for an arbitrary constant unit vector n. From this inequality we can conclude, that x-independent 
fluctuations of an exchange field of constant length around a ferromagnetic mean-field configuration 
do not decrease the free energy and hence are not a source of instability. 

Since the calculation is rather easy, we give an explicit formula for the coefficient Ct of the 
magnetic cr-model when the exchange field is ferromagnetically ordered. Because of the stability of 
the system with respect to x-independent fluctuations we must find a negative sign for this coefficient. 
We use the Holder inequality (|4.36|) . in order to get an upper bound for the large-/? behaviour of the 
trace in l)4.35j) and find, for large /3, and up to a W-independent constant 



S e ff(W) ~ {V- - V+) lndet(3 4 + W) + S K (W) 



(4.37) 
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where V± = Vol{p|e(p) + Wo/2 < 0}. In order to calculate Ct we take a look at the real part of 
lndet(d 4 + W), 

mndet {d t + W) = - bidet [(o» t + W)(-d t + W)] = - lndet (A + (d t W)) , 

where A := — <9 t 2 + (l/4)Wg . We treat (dfW) as a perturbation of A and proceed in the same way as 
in the last subsection to arrive at 

mndet (d t + W) = -^tr(A- 2 {d t W) 2 ) + h.o. , 

where we omitted W-independent terms. Finally, we find for the coefficient Ct 

w 2 g/_-y + ) v 1 (,w ^oo v--v+ 

* 8(2vr) 3 /? ^(kl + Wll^f 4(2vr) 3 ^ ' 1 ' 

which is indeed manifestly negative. 

Next, we investigate the stability of the ferromagnetically ordered mean-field configurations w.r.t. 
time-independent, but x-dependent perturbations, which amounts to the calculation of Cy in the 
magnetic u-model. This calculation is somewhat harder, because the p-dependent modes in the 
operator Z>w do not decouple. We perturb the ferromagnetic mean-field configuration W = Wqd. z 
by time-independent fluctuations, x> an d calculate their contribution to the real part of the effective 
action to second order in |x|/Wo. Since the effective action has a quadratic minimum w.r.t. the 
length of the exchange field, we can assume % t° De orthogonal to n z . Then (omitting independent 
terms) 



MS eff (W n z + X ) = -~ J> T (x)*r(x - y)x(y) - \ tr(A- x {e, x}^{e, X }) + h.o. , (4.39) 
where 



x.y 



9 9 W? 

A = -d 2 + e + + W ea 3 , 

In the trace we only keep contributions quadratic in |p|/Wo, where p is the momentum labelling 
the modes of x (derivative expansion). Spin degrees of freedom and Matsubara frequencies can be 
summed over explicitly, and - after a straightforward but tedious calculation - we obtain a formula 
for the coefficient, Cy, of the term 

U x 

in the expansion of the second term on the R.S. of Q4.39JI . which, for large values of /3, is given by 



4W 1 V 2^4 cosh/3(e(k)--^a) + l 



+ e('^-| 6 (k) 



PW ( 1 I 



v cosh (e(k) + Y) + 1 cosh ( e ( k ) - 2 ) + 1 , 
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By plugging W = Wbn z +Xi with x orthogonal to n z , into the a- model action (j4.21j) . one easily sees 
that the coefficient Cy just derived coincides with the corresponding coefficient in the cr-model. 

For the Hubbard model we can deduce sufficient conditions for the linear stability of the ferro- 
magnetic mean- field configuration (i.e., for Cy to be negative) from (|4.40f) . in the regime where (5 is 
large. We denote by Vp the Fermi velocities of the upper/lower spin band, namely 



Wq _i_ _i_ _|_ 

|e(k) ± -^-| ~ Vp\lt — k^l , for k close to the Fermi momenta k^, , 



by Ap the area of the Fermi surface of the upper/lower spin band, and by Vp the respective volumes. 
Then, for large /3, a sufficient condition for linear stability is that 

^(A+vp + A-Vp) > < ax (F- - V+) , (4.41) 

where v max = max{|Ve(k)| |e(k) - W /2 < < e(k) + W /2}. 

We now combine this stability condition with the mean-field equation (|4.7I) for the Hubbard- 
model, i.e., 

2U Vp - V+ 



to arrive at the stability condition 



U Ap:vp~ + ApVp 2 

> V 



3 (2vr) 3 

Details of these and other related calculations will appear in |18j . 



Magnon dynamics and the W Z-term 

In this subsection we consider the Wess-Zumino-term in Q4.21JI . We first consider an x-independent 
exchange field W(t) of constant length |W(i)| = Wq and calculate the variation of the first term on 
the R.S. of 1)3.40(1 for e = 0, given a variation, <5W, of the exchange field: 

(Hndetpo.w) = tr^wp^wA^wr^o.w) , (4-42) 

where 

5D ,w = ^5W(t) ■ a = 5W(t) . 

In order to evaluate the R.S. of (|4.42j) . we expand (Dq w Z?o,w) _1 m a Neumann series in powers of 
(d t W), 

(Dl^Do^r 1 = A- 1 + A-^WA" 1 + . . . , 
where A = — df + (1/A)Wq . Plugging this series into (|4.42|) . we find that 

5 In det(D ,w) = SS {1) (W) + 8S (2) (W) + . . . . (4.43) 

where 

5S (1) (W) = tridWA^D^) , (4.44) 

5S (2) (W) = tr {8W A' 1 {d t W) A' 1 D^ w ) , (4.45) 
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etc.. Using the cyclicity of the trace and the fact that Tr(<r) = (where "TV" is the trace on 2 x 2 
matrices), we see that 

P 



where 



5S^(W) = a(W ) J dtTr(W(t)5W(t)). 



and Zp = (tt//3)(2Z + 1). Since 

Tr(W(t)5W{t)) = ^5Tr(W(t) 2 ) , 

5S^ only depends on variations of the length of W(i) and vanishes if the length is held fixed. 
By tem 

5S (2) {W) = 5Sf\w) + 6S?}(W) , 

where 

5sf\w) = -tv(5WA~ 1 (d t W)A~ 1 d t ) = --^2tv(5W j A- 1 (d t W) j A- 1 d t ) 

3=1 

and 

SSffiW) = tr(5WA- 1 (5 t T^)A- 1 VF) = - % - J] e^trO^A" 1 ^)^ -1 ^) . 

j,l,m 

Here we have used that Tr (a j(Ti) = 25 ji and Tr(<7jcr;<7 m ) = —2ieji m . Since dt, A and Wj(t) are real 

(2) (2) 

operators, 5S} is real and 5Sj/ purely imaginary. The real part has been calculated in the last 

(2) 

subsection. Turning to 5Sj/ (W), we observe that 

5S$(W) = ti(W5W(d t W)A- 2 ) + higher derivative terms 



(Wo) / dt^e jlm W j (t)W(t)dtW m (t) + h.d.t. 



o 

I 



j,l,m 



o J dtW(t) ■ (5W(t) A d t W(t)) + O({f3W y l ) + h.d.t. , (4.46) 
o 

where, in the last equation, we have used that 

We note that h.d.t. ~ O^Wq)- 1 ). 

Recall that W(i) has periodic boundary conditions at t = 0, /3. Viewing the imaginary-time circle 
[0, j3) as the equator of a sphere of radius (5/ (2tt), we may extend the unit-vector field W(i) from 
the equator of 5J to a continuous unit-vector field W(i,s), < s < (3/2, on the entire sphere in an 
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arbitrary way, but with W(i,/3/4) = W(i). The first term on the R.S. of ()4.46|) then turns out to 
be the variation of the functional 



0/4 



Swz(W) = -^ J dt j dsW(t,s)-(d s W(t,s)Ad t W(t,s)) 
o 

W • (dW A dW) , (4.48) 



o o 

i 



2 

5 2,JV 

where "d" denotes the exterior derivative, and S 2 ' N is the northern hemisphere of S 2 . We recognize 
(|4.48|) to be the Wess-Zumino term in (j4.21j) for an x-independent (i.e., ferromagnetically ordered) 
exchange field. The corresponding coefficient CVz depends on e and will be determined later on. 
The R.S. of ()4.48j) is expressed in standard differential-form notation, which makes it manifest that 
Swz(W) is independent of the radius, 0/2n, of the sphere. 

The leading term on the R.S. of (|4.46|) is also the variation of 



S'wz(W) = \ f W • (dW A dW) , 



52,s 

where S 2 ' S is the southern hemisphere of S 2 . One observes that 

Swz(W) -S' wz (W) = -^Jw- (dW A dW) = 2irin , (4.49) 

s 2 

where n = n(W) £ Z is the degree of the map W : S 2 — > S 2 . It follows that, apparently, the 
Wess-Zumino action Sy/Z is only determined modulo an integer multiple of 2iri. Due to l|4.49|) ■ 
exp(5pvz(W)) is a single-valued functional of W. 

If the exchange field W(x,t) is antiferromagnetically ordered, i.e., the Neel field 

N(x,t) = (-l) |x| W(x,t) = W n 



is a constant field pointing in the direction of some unit vector n, the Wess-Zumino term in (|4.21|) 
obviously tends to in the formal continuum limit. It follows that, in antiferromagnetically ordered 
states, at sufficiently low temperatures, there are gapless Goldstone bosons with dispersion 

w(k)wi;|k|, (|k|«0), (4.50) 

for some velocity v. 

In ferromagnetically ordered systems, however, the Wess-Zumino term does not disappear. Using 
(|PT)) and (|4~15|) . the coefficient C W z in ijOTj) is found to be 

°™ = -ri&r- (4 - 51> 

in the limit where (3Wq — > oo. Had we normalized the Fourier modes and were we in a finite volume, 
we would find that exp(S* e jj(W)) is single- valued. From (|4.21l) we derive the equations of motion 

C wz W x A^V x = J2K(x-y)W y + aW x , (4.52) 
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with ~W x (t) = W(x, t), W x = dtW x , where K(x — y) = WqK(x — y) + C-sj8\ x -y\\, and where a is a 
Lagrange multiplier arising from the constraint |W x (t)| 2 = 1, for all x and all t. Taking the vector 
product of (|4.52j) with W x , and doing the Wick rotation from imaginary time back to real time |23) . 
we obtain the equations of motion for magnons in a ferromagnet 

iCwz&x =W X AJ2 K(x - y)W y , (4.53) 

y 

which is the well known Landau- Lifshitz equation] see also |19j . 

If J, and hence K, are of short range and if the equilibrium state of the system were ferromag- 
netically ordered, i.e., 

W(x,t) =W n + X (x,t), W >0, 

where n is a unit vector and x( x i t) is assumed to be small, then the dispersion of magnons is found 
to be 

Ikl 2 

wfk) = — — , for some constant M > , 
as expected. If n points in the z-direction then x hes essentially in the x — y plane. We define 

<t>(p,t) = Xi{x,t) +ixi(x,t) . 

Passing to the formal continuum limit, the action functional for the complex magnon field <j) of a 
ferromagnet is seen to be given by 

5 mag „on(0,0) = J dtj dx$(x,t)d t (f)(x,t) + ^4>(x,t)(Acj))(x,t)] , (4.54) 



which is the action for a system of conserved non-relativistic bosons with vanishing chemical potential. 
Since W\ = Xi = (0+ l ^)/2 is an observable field, the chemical potential actually necessarily vanishes, 
and there is no "magnon condensation"; (just like photons do not form any Bose condensates). 

Absence of symmetry breaking in one and two dimensions at positive tempera- 
tures 

For (3 < oo, we can study the fluctuations of the modes W(&o,k) of the exchange field with ko = 0, 

P 

i.e., of JW(x,t)dt. 
o 

We set 

P 

J W(x,t)dt = l3(W m f + xi(x)n z + X t(x)), 
o 

where W" 1 -^ is a solution of the ferromagnetic or the antiferromagnetic mean-field equation. Then, 
in d spatial dimensions, we conclude from our stability analysis 

(Xt(x) 2 ) f3 4l OC 

which, for d = 1 and 2, is infrared-divergent. Thus, transversal fluctuations around the mean field 
are gigantic and hence destroy long-range order; (this is the Mermin- Wagner theorem, see, e.g., 0). 
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In two dimensions, there is an alternative way to understand the absence of ordering when 
P < oo: Setting W(x,t) = W(x) (ferromagnetic short-range order), or W(x,i) = (— l)^ W(:r) 
(antiferromagnetic short-range order), where W(x) is slowly varying, and passing to the formal 
continuum limit (lattice spacing — > 0), one finds that the action 5 e //(W) has (approximate) critical 
points, indeed local maxima, on time-independent configurations W(x) with the properties: 

|W(s)|wW , for all x, 

where Wq is a solution of the (ferromagnetic or antiferromagnetic) mean-field equation, 

W(x) — > Won , as |x| — > oo , 

where n is an arbitrary unit vector, and the degree of the map W(x), as measured by the integer 
("winding number") 

1 f j2 «r/ \ ( 9W(x) dW{x)\ 

n / d 2 xW(x) ■ n v ' A - v ; , 

4ttW$ J V ' \ dxi dx 2 J 

is non-zero. Such configurations are called "instantons" . The contributions of instanton configura- 
tions to a functional integral, such as (|3.43|) . destroy long-range order. 

Remark. Generally speaking, it is hard to justify the use of steepest descent in approximately 
evaluating functional integrals, such as (|3.41j) and l|3.43|) . because there isn't any large constant 
N = h~ 1 multiplying S e ff(W). If, however, we place N identical species of spin-^-fermions on 
each site x £ T, coupled to each other only through exchange interactions between the total spin 
operators, and if we set Kq = MKq then S e ff(W) = S^^(W) is replaced by AfS^\"W), and 
steepest descent becomes reliable. 



5 Reflection positivity and phase transitions in Heisenberg-models 

In order to establish the existence of phase transitions in some of the models introduced in Section 
3, we recall the method of reflection positivity (|5|.|l].|5n|.[21|). F° r simplicity, we choose the lattice 
T to be given by 7L d . We decompose T into two equal halves 

r + = {x g r\xi > o}, r_ = {x e r|xi < o}, (5.1) 

where Xj is the j th component of x G T. Let II be the plane in R orthogonal to the 1-direction at 
height xx = h, and let i? denote reflection at II. Clearly, i? maps T_ onto r+, and conversely. 
An exchange coupling matrix, J(x,y), is said to be reflection positive iff 

J2J(x)J(x,'&y)f(y)>0, (5.2) 

x,y 

for all functions / on T with supp / C T + . Note that this condition is independent of the diagonal 
elements, J(x,x), of J; so, by choosing them appropriately, we can always ensure that J is either 
positive-definite or negative-definite. By computing some integrals with a Gaussian measure of mean 
and covariance — J -1 (J negative-definite) we see that if J is reflection positive then so is — J -1 ; 
see e.g. |2U|.|21|. Furthermore, if J is reflection-positive then so is — J, where J is given by 

J(x,y) = (-l)\*+y\ J(x,y) = (-l)W+Mj(x,y) (5.3) 
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with \x\ = Sj=i x j- (Note that if J is positive-definite then J is positive-definite, too.) The following 
simple calculation proves our claim. Let 

/» = (-l)W/(x). 
Since J is assumed to be reflection-positive, 

< ^J(x)J(x,tiy)f(y) 

= f 7R(-i) w+l ^(^^)/(2/) 

x,y 

= -^7M(-l)I^J(x,^)/(y) 

x ' y _ 

x,y 

for an arbitrary function / with supp/ C T + (hence supp/ C These considerations are 

summarized as follows: 

If J is reflection-positive then —K := —J^ 1 is reflection-positive, too; (5-4) 

and 

if— J is reflection-positive (i.e., J is reflection-negative) 

then —K := — J -1 is reflection-positive, too. (5-5) 

In both situations, we can choose the diagonal elements of J such that J and K are positive-, or 
negative-definite. 

A reflection-positive matrix J is called ferromagnetic, while a reflection-negative matrix is called 
antiferromagnetic. It is easy to see that, for nearest-neighbour exchange couplings, these notions of 
"ferro-" and "antiferromagnetic" coincide with the familiar ones. 

Let J be ferromagnetic and negative-definite. Then the Gaussian measure 

dn-jQ) = const ■e^*,y* {x)J{x > v) * {y) Y[d n <t>(x) : 

X 

where <f>(x) 6 M n , n = 1,2,3,..., and the constant is chosen such that d\i_j is normalized, is 
reflection-positive, in the sense that 

J dn-j(4>)(QF)&)F&) > 0, (5.6) 

for an arbitrary bounded function F that depends only on the variables {<j>(x)\x S r + }; the space of 
all such functions is denoted by J-+. The function QF in (|5.6[) is defined by 

(GF)(<t>)=F(fo), (5.7) 

where <f>&(x) = </>($x). If <f>^ , ■ ■ ■ ,<f>^ are independent Gaussian random fields with distributions 
given by dfi_ J (t)(^ t '), where is ferromagnetic and JW < 0, t = 1, . . . , I, then the distribution 

l 

d^)=II^-rWjw(^ (t) ) ( 5 - 8 ) 
t=l 
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of the random field (f> given by 

<t>(x,t)=<f>V(x),t = l,...,l, 

is reflection-positive, for arbitrary > 0, t = 1, as is easily seen. 

These considerations have the following consequence that will play a crucial role in our analysis of 
phase transitions: Let W(x,t), x E T, t E [0, (3), with W(x,t + /3) = W(x,t), be the imaginary time 
exchange field, and let d//#-(W) = exp(S'^:(W))DW be its Gaussian distribution, as introduced 
in Section 03 1)3. 2 7 J) through 1)3.31)1 . If the exchange couplings J are ferromagnetic (in the sense 
specified above) and positive-definite then —K = — J -1 is ferromagnetic and negative-definite; hence 
the Gaussian measure 

d[iK(W) is reflection-positive, (5-9) 

as follows from (f3~2"9"|) . JOJ and (Oil . 

Next, we suppose that the exchange couplings J are antiferromagnetic (in the sense specified 
above). Then, by 1)5.5)1 . —If = —J" 1 is reflection-positive and negative-definite, for J positive- 
definite. We introduce a random field, N, ("N" for "Neel"), by setting 

N(x,t) = (-l) |x,| W(x,t). (5.10) 

If W has distribution dfiK (W), with if = J^ 1 , then N has distribution d//g.(N), and hence 

dfi^.(N) is reflection-positive. (5-H) 

We set 

Ntf(x,t) = N(tfcc,t), (5.12) 

and (see 1)577)1 ) 

(9F)(N) =F(N„), (5.13) 

for an arbitrary function F of {N(x, i)|a: E T, t E [0, /?)}. 

We consider the models introduced in Section |21 without hopping term but with M Fermion 
species. They have an action functional given by 



S(C, C; W) = ][; / dtJ2 \c { s a) (x,t)(^-A C< a > (x,t)- u a W(x, t) • S ^ (x, t) 

a=l q x L V / 

-\\ dt(j2 w ( x ^) K ( x -y) w (y^)] > ( 5 - 14 ) 

\x,V / 



where v a = ±1, for all a = 1, . . . and K(k) = Uu + J^{k)j ; see ()3~20)l . ()3~29l . and (EQ2l) 

through 1)3.34)1 . Since the first term on the R.S. of 1)5.14)1 does not couple different sites in the lattice, 
the effective action of the exchange field W, see ()3.36)) and ()3.4U|I . is given by 



(W) = -^ I dtlj2W(x,t)K(x-y)W(y,t)\ + ££lndet(A, aW(a , ) 

V x,y / ax 



(5.15) 
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with Z?w given in (13.381) . for t(x) = 0. We define 



un) = n det (^(-DWN(x,^ 



a=l 

N 



(5.16) 



Y[ det(D VoW(!B) .)), 



0=1 



Then it follows that 



e/ x (N) = /x(N^(x,-)) 

= n det (^(-i)i-iN^ 



a=l 
A/" 



n det ( D -".(-i)i**iN(*x,o) ( 5 - i? ) 



a=l 

J\T 



n det (^ a 



(-l)l<MN(0ay),' 



a=l 



/<to(N), 



where, in the third equation, we have used that (— 1)^ = — (— l)'* 1 ', and, in the fourth equation, we 
have inserted identity Q4.27JI . Equation (|5. 17|) implies that 



n /*( n )) = n /^ n )- ( 5 - i§ ) 



xer + x£T- 

From this equation and equations (|5.11|) and Q5.6JI we conclude the following simple, but quite fun- 
damental 

Result A Suppose the exchange couplings J are antiferromagnetic (in the sense introduced above) 
and positive-definite. We set 

d Mf (N) = const • Yl f x (N)d^(N) (5.19) 

X 

= const' • exp (S e// (N))pN, (5.20) 
where the constants are chosen such that J d(i e ff(N) = 1, and 

Se//(N) ■= ~\J dt(TN(x,t)K(x-y)N(y,t)\ 



(5.21) 



lN(av) 



with K = J" 1 ; see (t5~T5l . (t5~TTH) . (t5~Tll . 

Then the measure dfi e ff(N) is reflection-positive, in the sense that 

y d/i e// (N)(6F)(N)F(N) > 0, (5.22) 
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for arbitrary functions F € J-+ (where J-+ is denned right above (|5.7|) ). 



Next, we suppose that the number of fermion species, N = 2A4 is even, with 

ui = ... = v M = 1, VM+i = ■■■ = v 2 M = -1- (5-23) 

We define 

g x (W) = det(D WM ) M det(D_ WM ) M . 
Identity (|4.27|) then implies that g x (W) > 0, for all x £ T, and, setting 

W#(x,t) =W(dx,t), 

we find that 

e^(W) = @g x ( W) = g# x ( W) . (5.24) 

Defining 

d/x e //(W) = const- n^(W)^(W) (5.25) 

a; 

= const' -exp (%(W))pW, (5.26) 

where K = J -1 , S e ff(W) is as in (|5.15|) . and the constants are chosen such that J dfi e ff(W) = 1, 
we arrive at 

Result F Suppose the exchange couplings J are ferromagnetic (in the sense introduced above). 
Then the measure dfj, e ff(W) is reflection-positive, i.e. 

J d f i eff (W)(QF)(W)F(W) > 0, (5.27) 

for arbitrary functions F £ T+. 

From Results A and F one obtains the following infrared (spin-wave) bounds, originally discov- 
ered in |Sj and generalized in [2"U].|21|. 

(IRA ) Under the hypotheses of Result A (in particular, J antiferromagnetic), one has the inequalities 

< (N(k,fc )-N(-k,-fc )W 

/- ~ \-i (5.28) 



for some constant Wq > 0. Note that N(— k, — ko) = N(k, ko), because N(x,t) is real. 

If Wq is strictly positive then the state ((•)) ( a j ^ exhibits long-range order, and if we couple N to 
an arbitrarily small external field, en, |n| = 1, then 

lim(N(x, i)W„ = W n. 3 (5.29) 



3 This claim is heuristic, but can be replaced by an equivalent, mathematically rigorous one, |20|. 
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Inequality (|5.28|) implies that, for an arbitrary bounded function f(x, t) of rapid decay in x £ T and 
periodic in t £ [0, fi), 



< (|N(/)| 2 ) AAt 



< /W 2 |/(0,0)| 2 + £ I 



dkJ /(k '**T , (5-30) 



2 



fcoe - z ^ r K(k)-K(0) 



where N(/) = / dt £ N(a, i)/(a;, t). 

as ___ 

Returning to identities (|3.42|) and (|3.43[l . with 



A/" 

S(x,t) = Y,^ a \x,t), (5.31) 



a=l 

and recalling that W(x,i) = (— l)l :r lN(x, t), we conclude from (j5.3Uf) that 



0<((fs(0,t)dt) 2 )^ (5.32) 



o 

\2 



< pK(0) + /3W 2 K(0) 2 + I dk ^ K( ^\ — . (5.33) 

I K(k) - K(0) 

We recall that —K is reflection-positive, and the diagonal elements are chosen such that K (or, 

equivalently, K and J) is positive-definite. Therefore K(k) > 0, for all k G B^. Furthermore, 
if J is of short range then K is of short range, too, i.e., K(x) decreases rapidly in \x\, and, by 

reflection-positivity of —K, K(0) < K(k). It follows that 

< K(k) -K(0) < K|k| 2 , (5.34) 

for some constant k. Thus, under the assumptions just specified, the third term on the R.H.S. of 
(|5.32j) is bounded above by a finite constant, I{K, d), in dimension d > 3, but is (infrared-) divergent 
in one and two dimensions; (I(K,d = 1,2) is finite only if J is of long-range; see |2*0].|21|). 

A fairly elementary inequality due to Bruch and Falk, see [1], yields the following lower bound 

P 

<( J S(0,t)dt) 2 )^ > C((3^,M)f3 (5.35) 
o 

for some constant C((3, (a,N) > 0. If the on-site repulsion Uq, the chemical potential \x and the value 

of if _1 (0) = f B dkK -1 ^) are tuned appropriately such that, at zero temperature, the band of each 
species of fermions is half-filled then C((3,[i,N) > and limjv"->oo C(/3, (J,, A/) = oo. 

Thus, we conclude that, for appropriate choices of Uq, (jl and K (0), and for a sufficiently large 
number, N , of fermion species, 

Wq is strictly positive, for (3 large enough (5.36) 
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Hence, by (I3~i2l and 



lim(S(x,0))^ i£n = -^K(x-y)lim(W(y,0)) A/ , j£n 

y 

= -(-l) N ^^-y)lim(N(y,0)) AMi£n (5.37) 

= -(-l) N K(0)TU n / 0, 

with M = K(0)\Wq\ the local magnetization, i.e., the system exhibits Neel order, for sufficiently 
large values of (3 and J\f. A similar result has first been proven in 0]. 

Arguments related to those used in the proof of the Mermin- Wagner theorem (see jH]) show that 
if M > then there are gapless spin waves (see e.g. |24j ) . in accordance with the mean-field picture 
of Section 0J 

Next, we exploit Result F. It implies that, under condition ()5.23j) . the following infrared bounds 
hold. 



(IRF) Under the hypotheses of Result F (in particular, J ferromagnetic and ()5.23|) holds), the 
inequalities 

< <W(k,fc )-W(-k,-&o)W 

-i (5.38) 



< pW Q %(k)5 ko>0 + (k) - K(0) 

for a constant IU 2 > 0, are valid. (See , H3 , EU ) • The consequences of these inequalities are 
perfectly analogous to those discussed for antiferromagnets above. We set 

M 

S(x,t) = Y, (s (a) M) - S( a+M \x,t)), (5.39) 

a=l 

see ()5.23|) . Then, under appropriate conditions on Uq, /x and K _1 (0) = J(0), and for d > 3, 

M = lim(S(x, 0)^, £ n = -£(0) W n ± (5.40) 

for (3 and A4 large enough, which proves long-range order at low temperatures, for M large enough; 
(presumably N = 2M. = 2 suffices if J is nearest-neighbour and d > 3). 

Unfortunately, there is, as yet, no such result on the existence of a phase transition in a quantum 
Heisenberg ferromagnet with only one species of fermions (electrons), i.e., J\f = 1. The reason is that 
this model is not described by a reflection-positive distribution, dfi e ff(W), for the exchange field 
W. Setting 

f x (W) = det(D WM ), (5.41) 

dfj, e ff is given by 

d/i e //(W) = const- J]/,(W)d^(W). (5.42) 

a; 

Now, if J is ferromagnetic d^i^(W) is reflection-positive if we set W#(x,i) = W(i?x,t). However, 
using (|4.27j) . 

©/*(W) = det(Z) w( ^ i . ) ) 
= det(Z?_ W (^ v) ) 

= /*«(-w) ^ U X (W). 
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Thus, dfx e ff(W), as given in (|5.42[) . is not reflection-positive, and there is therefore no reason why 
the Infrared Bounds (IRF) should hold. A very similar problem is encountered in the study of 
Bose-Einstein condensation in lattice models of non-relativistic, interacting bosons. 

In order to get some preliminary insights into phase transitions for the Heisenberg ferromagnet, 
one may consider the following "static approximation": one replaces /^(W) by 

/J°)(W) = det(I} wC o) (a;) ), (5.43) 

where 

i r 13 

W (0) (z) = - / dtW(x,t) 
P Jo 

which is time-independent. Then 



e/r (W) = /r(-W) = /^(W), (5.44) 
as is easily checked, and /i (W) > 0. Replacing d/i e ff(W) by 

^° / ) / (w)=n/i° ) (w)^(w), 

X 

we conclude that d//^j(W) is reflection-positive if J is ferromagnetic, hence the infrared bounds 
(|5.38f) hold, and we conclude that, in the static approximation, for appropriate choices of Uq, /j, and 
J(0), a phase transition accompanied by continuous symmetry breaking occurs at sufficiently low 
temperatures. 

If one replaces electrons with spin ^ by fermions with spin s, or by 2s identical species of spin-i 
fermions (with u a = 1, for all a = l,...,2s), and if one then takes s — ► oo (rescaling the spin 
operators, S(x,t) iS(x,i)) then the static approximation becomes exact. This is not surprising, 
because the limit s — > oo corresponds to the classical limit, as shown in [25], and the spin-s Heisen- 
berg model approaches the classical Heisenberg model for which the existence of a phase transition 
has been established in [SJ. Our formalism, involving the exchange field W, offers a neat and simple 
way of recovering some of the results in [23]; (we leave this as an exercise to the reader). 

Remarks. 

1. We have outlined a proof of existence of a phase transition for a class of Heisenberg antiferro- 
magnets; (the proof is mathematically rigorous; some missing details can be inferred from [3], 
0]). Of course, this is not a surprise. In [I], Dyson, Lieb and Simon have already proven a 
similar result also using reflection positivity (albeit in a different manner). Their method of 
proof, too, breaks down for quantum ferromagnets. 

2. It appears that to understand phase transitions in quantum ferromagnets and the related 
phenomenon of Bose-Einstein condensation for lattice gases of non-relativistic bosons in a 
mathematically rigorous way, we would have to resort to a full-fledged renormalization group 
analysis - a technically rather demanding task. 

3. The fact that the function f x (W) defined in (|5,41|) (see also (|5.19[) ) is neither positive, nor even 
real is the origin of the "sign (complex phase) problem" in numerical simulations of quantum 
ferro- and antiferromagnets involving, e.g., the Monte Carlo method applied to (|5.42[) . The 
great quality of methods based on reflection positivity is that they are applicable to the analysis 
of certain complex distributions ^Me//- 



= iw(x,fc Q = 0)), 
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4. If one considers models in d dimensions with SO(2n) spins, n = 1,2,3,..., then for d > 3 one 
can prove a ferromagnetic phase transition, using arguments very similar to those described 
above. 

5. After completion of this paper, work of Bach et al. [57] appeared to which we draw the reader's 
attention. 
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